Drivers of soil microbial and detritivore activity across global grasslands.

Covering approximately 40% of land surfaces, grasslands provide critical ecosystem services that rely on soil organisms. However, the global determinants of soil biodiversity and functioning remain underexplored. In this study, we investigate the drivers of soil microbial and detritivore activity in grasslands across a wide range of climatic conditions on five continents. We apply standardized treatments of nutrient addition and herbivore reduction, allowing us to disentangle the regional and local drivers of soil organism activity. We use structural equation modeling to assess the direct and indirect effects of local and regional drivers on soil biological activities. Microbial and detritivore activities are positively correlated across global grasslands. These correlations are shaped more by global climatic factors than by local treatments, with annual precipitation and soil water content explaining the majority of the variation. Nutrient addition tends to reduce microbial activity by enhancing plant growth, while herbivore reduction typically increases microbial and detritivore activity through increased soil moisture. Our findings emphasize soil moisture as a key driver of soil biological activity, highlighting the potential impacts of climate change, altered grazing pressure, and eutrophication on nutrient cycling and decomposition within grassland ecosystems.

Multidimensional responses of grassland stability to eutrophication.

Eutrophication usually impacts grassland biodiversity, community composition, and biomass production, but its impact on the stability of these community aspects is unclear. One challenge is that stability has many facets that can be tightly correlated (low dimensionality) or highly disparate (high dimensionality). Using standardized experiments in 55 grassland sites from a globally distributed experiment (NutNet), we quantify the effects of nutrient addition on five facets of stability (temporal invariability, resistance during dry and wet growing seasons, recovery after dry and wet growing seasons), measured on three community aspects (aboveground biomass, community composition, and species richness). Nutrient addition reduces the temporal invariability and resistance of species richness and community composition during dry and wet growing seasons, but does not affect those of biomass. Different stability measures are largely uncorrelated under both ambient and eutrophic conditions, indicating consistently high dimensionality. Harnessing the dimensionality of ecological stability provides insights for predicting grassland responses to global environmental change.

Herbivory and nutrients shape grassland soil seed banks.

Anthropogenic nutrient enrichment and shifts in herbivory can lead to dramatic changes in the composition and diversity of aboveground plant communities. In turn, this can alter seed banks in the soil, which are cryptic reservoirs of plant diversity. Here, we use data from seven Nutrient Network grassland sites on four continents, encompassing a range of climatic and environmental conditions, to test the joint effects of fertilization and aboveground mammalian herbivory on seed banks and on the similarity between aboveground plant communities and seed banks. We find that fertilization decreases plant species richness and diversity in seed banks, and homogenizes composition between aboveground and seed bank communities. Fertilization increases seed bank abundance especially in the presence of herbivores, while this effect is smaller in the absence of herbivores. Our findings highlight that nutrient enrichment can weaken a diversity maintaining mechanism in grasslands, and that herbivory needs to be considered when assessing nutrient enrichment effects on seed bank abundance.

Clarifying the effect of biodiversity on productivity in natural ecosystems with longitudinal data and methods for causal inference.

Causal effects of biodiversity on ecosystem functions can be estimated using experimental or observational designs - designs that pose a tradeoff between drawing credible causal inferences from correlations and drawing generalizable inferences. Here, we develop a design that reduces this tradeoff and revisits the question of how plant species diversity affects productivity. Our design leverages longitudinal data from 43 grasslands in 11 countries and approaches borrowed from fields outside of ecology to draw causal inferences from observational data. Contrary to many prior studies, we estimate that increases in plot-level species richness caused productivity to decline: a 10% increase in richness decreased productivity by 2.4%, 95% CI [-4.1, -0.74]. This contradiction stems from two sources. First, prior observational studies incompletely control for confounding factors. Second, most experiments plant fewer rare and non-native species than exist in nature. Although increases in native, dominant species increased productivity, increases in rare and non-native species decreased productivity, making the average effect negative in our study. By reducing the tradeoff between experimental and observational designs, our study demonstrates how observational studies can complement prior ecological experiments and inform future ones.

Environmental heterogeneity modulates the effect of plant diversity on the spatial variability of grassland biomass.

Plant productivity varies due to environmental heterogeneity, and theory suggests that plant diversity can reduce this variation. While there is strong evidence of diversity effects on temporal variability of productivity, whether this mechanism extends to variability across space remains elusive. Here we determine the relationship between plant diversity and spatial variability of productivity in 83 grasslands, and quantify the effect of experimentally increased spatial heterogeneity in environmental conditions on this relationship. We found that communities with higher plant species richness (alpha and gamma diversity) have lower spatial variability of productivity as reduced abundance of some species can be compensated for by increased abundance of other species. In contrast, high species dissimilarity among local communities (beta diversity) is positively associated with spatial variability of productivity, suggesting that changes in species composition can scale up to affect productivity. Experimentally increased spatial environmental heterogeneity weakens the effect of plant alpha and gamma diversity, and reveals that beta diversity can simultaneously decrease and increase spatial variability of productivity. Our findings unveil the generality of the diversity-stability theory across space, and suggest that reduced local diversity and biotic homogenization can affect the spatial reliability of key ecosystem functions.

Successional patterns of terrestrial wildlife following deglaciation.

Disturbance is a key driver of community assembly and patterns of diversity. Whereas successional changes in vegetation have been well-studied, postdisturbance successional patterns of wildlife communities remain poorly understood. Here, we investigated the roles of site age and habitat in shaping community assembly and the diversity of terrestrial mammals in Glacier Bay National Park, Alaska (GBNP), which has undergone the most rapid and extensive deglaciation in the world since the Little Ice Age. Deglaciation has extensively altered the landscape, opening up new habitat for recolonization by plants and animals. We used camera traps, small mammal trapping and vegetation surveys to investigate the patterns of mammalian succession and beta diversity following deglaciation, using a space-for-time substitution across 10 sites during summers 2017 and 2018. Site age and habitat characteristics were not strongly correlated (r < 0.46), allowing the influences of time since disturbance and habitat changes to be distinguished. PERMANOVA analyses indicated that mammal community assembly was more strongly influenced by site age than habitat, whereas habitat and age had similar effects on beta (between site) diversity. Beta diversity was higher for smaller, less mobile mammals than larger, more mobile mammals and was primarily driven by species turnover among sites, whereas relative turnover was much lower for larger mammals. A comprehensive review of historical distributions of mammals in GBNP supported our findings that species turnover is a driving influence of community assembly for smaller mammals. Our results indicate that body size of mammals may play an important role in shaping colonization patterns postdisturbance, likely via size-related differences in mobility. Patterns of wildlife community assembly may therefore not track vegetation succession following disturbances if there are barriers to movement or if dispersal ability is limited, highlighting the importance of incorporating landscape connectivity and species traits into wildlife conservation efforts following disturbances. This knowledge may improve predictions of mammalian community assembly following major disturbance events.

Linking changes in species composition and biomass in a globally distributed grassland experiment.

Global change drivers, such as anthropogenic nutrient inputs, are increasing globally. Nutrient deposition simultaneously alters plant biodiversity, species composition and ecosystem processes like aboveground biomass production. These changes are underpinned by species extinction, colonisation and shifting relative abundance. Here, we use the Price equation to quantify and link the contributions of species that are lost, gained or that persist to change in aboveground biomass in 59 experimental grassland sites. Under ambient (control) conditions, compositional and biomass turnover was high, and losses (i.e. local extinctions) were balanced by gains (i.e. colonisation). Under fertilisation, the decline in species richness resulted from increased species loss and decreases in species gained. Biomass increase under fertilisation resulted mostly from species that persist and to a lesser extent from species gained. Drivers of ecological change can interact relatively independently with diversity, composition and ecosystem processes and functions such as aboveground biomass due to the individual contributions of species lost, gained or persisting.

Root hemiparasitic plants are associated with more even communities across North America.

Root hemiparasitic plants both compete with and extract resources from host plants. By reducing the abundance of dominant plants and releasing subordinates from competitive exclusion, they can have an outsized impact on plant communities. Most research on the ecological role of hemiparasites is manipulative and focuses on a small number of hemiparasitic taxa. Here, we ask whether patterns in natural plant communities match the expectation that hemiparasites affect the structure of plant communities. Our data were collected on 129 national park units spanning the continental United States. The most common hemiparasite genera were Pedicularis, Castilleja, Krameria, and Comandra. We used null models and linear mixed models to determine whether hemiparasites were associated with changes in community richness and evenness. Hemiparasite presence did not affect community metrics. Hemiparasite abundance was positively associated with increasing evenness of herbaceous species, but not with species richness. The associations that we observed on a continental scale are consistent with evidence that the impacts of root hemiparasitic plants on evenness can be substantial and abundance dependent but that effects on richness are less pronounced. Hemiparasites mediate competitive exclusion in communities to facilitate species coexistence and merit consideration of inclusion in ecological theories of coexistence.

Evolutionary history of grazing and resources determine herbivore exclusion effects on plant diversity.

Ecological models predict that the effects of mammalian herbivore exclusion on plant diversity depend on resource availability and plant exposure to ungulate grazing over evolutionary time. Using an experiment replicated in 57 grasslands on six continents, with contrasting evolutionary history of grazing, we tested how resources (mean annual precipitation and soil nutrients) determine herbivore exclusion effects on plant diversity, richness and evenness. Here we show that at sites with a long history of ungulate grazing, herbivore exclusion reduced plant diversity by reducing both richness and evenness and the responses of richness and diversity to herbivore exclusion decreased with mean annual precipitation. At sites with a short history of grazing, the effects of herbivore exclusion were not related to precipitation but differed for native and exotic plant richness. Thus, plant species' evolutionary history of grazing continues to shape the response of the world's grasslands to changing mammalian herbivory.

Cover and density of southwestern ponderosa pine understory plants in permanent chart quadrats (2002-2020).

This data set consists of 101 permanent 1 m × 1 m (1-m2 ) quadrats located within southwestern ponderosa pine ecosystems near Flagstaff, Arizona, USA. Individual plants in these quadrats were identified and mapped annually for 19 years (2002-2020). The original chart quadrats were established between 1912 and 1927 to determine the effects of domestic livestock grazing on herbaceous plants and pine seedlings. Today these data provide opportunities to examine the effects of climate and land-use change on plant demography, population dynamics, and community processes. We provide the following data and data formats: (1) digitized maps of all plant locations in shapefile and geodatabase format, (2) shapefiles showing annual locations of each individual plant species, (3) annual maps of each quadrat in TIFF and PDF format, (4) annual basal area of each species per quadrat for species mapped as polygons, (5) tabular representation of polygon areas and centroid locations for plant species mapped as polygons, (6) tabular representation of point locations for plant species mapped as points, (7) plot-scale 20 m × 20 m overstory tree canopy cover, tree basal area, parent material, and elevation, (8) quadrat-scale information (e.g., quadrat site and number, coordinates in UTM Zone 12 and latitude/longitude, and descriptive comments for each quadrat), (9) plant species list, (10) summary of plant species observed in each quadrat, (11) summary of quadrats mapped by site and year, and (12) data formatted for use in Integral Projection Models (IPM) and plant population analyses. There are no copyright restrictions; please cite this paper and the associated data set when data are used in publications.

Nutrients and herbivores impact grassland stability across spatial scales through different pathways.

Nutrients and herbivores are well-known drivers of grassland diversity and stability in local communities. However, whether they interact to impact the stability of aboveground biomass and whether these effects depend on spatial scales remain unknown. It is also unclear whether nutrients and herbivores impact stability via different facets of plant diversity including species richness, evenness, and changes in community composition through time and space. We used a replicated experiment adding nutrients and excluding herbivores for 5 years in 34 global grasslands to explore these questions. We found that both nutrient addition and herbivore exclusion alone reduced stability at the larger spatial scale (aggregated local communities; gamma stability), but through different pathways. Nutrient addition reduced gamma stability primarily by increasing changes in local community composition over time, which was mainly driven by species replacement. Herbivore exclusion reduced gamma stability primarily by decreasing asynchronous dynamics among local communities (spatial asynchrony). Their interaction weakly increased gamma stability by increasing spatial asynchrony. Our findings indicate that disentangling the processes operating at different spatial scales may improve conservation and management aiming at maintaining the ability of ecosystems to reliably provide functions and services for humanity.

Nutrient identity modifies the destabilising effects of eutrophication in grasslands.

Nutrient enrichment can simultaneously increase and destabilise plant biomass production, with co-limitation by multiple nutrients potentially intensifying these effects. Here, we test how factorial additions of nitrogen (N), phosphorus (P) and potassium with essential nutrients (K+) affect the stability (mean/standard deviation) of aboveground biomass in 34 grasslands over 7 years. Destabilisation with fertilisation was prevalent but was driven by single nutrients, not synergistic nutrient interactions. On average, N-based treatments increased mean biomass production by 21-51% but increased its standard deviation by 40-68% and so consistently reduced stability. Adding P increased interannual variability and reduced stability without altering mean biomass, while K+ had no general effects. Declines in stability were largest in the most nutrient-limited grasslands, or where nutrients reduced species richness or intensified species synchrony. We show that nutrients can differentially impact the stability of biomass production, with N and P in particular disproportionately increasing its interannual variability.

Temporal rarity is a better predictor of local extinction risk than spatial rarity.

Spatial rarity is often used to predict extinction risk, but rarity can also occur temporally. Perhaps more relevant in the context of global change is whether a species is core to a community (persistent) or transient (intermittently present), with transient species often susceptible to human activities that reduce niche space. Using 5-12 yr of data on 1,447 plant species from 49 grasslands on five continents, we show that local abundance and species persistence under ambient conditions are both effective predictors of local extinction risk following experimental exclusion of grazers or addition of nutrients; persistence was a more powerful predictor than local abundance. While perturbations increased the risk of exclusion for low persistence and abundance species, transient but abundant species were also highly likely to be excluded from a perturbed plot relative to ambient conditions. Moreover, low persistence and low abundance species that were not excluded from perturbed plots tended to have a modest increase in abundance following perturbance. Last, even core species with high abundances had large decreases in persistence and increased losses in perturbed plots, threatening the long-term stability of these grasslands. Our results demonstrate that expanding the concept of rarity to include temporal dynamics, in addition to local abundance, more effectively predicts extinction risk in response to environmental change than either rarity axis predicts alone.

Species loss due to nutrient addition increases with spatial scale in global grasslands.

The effects of altered nutrient supplies and herbivore density on species diversity vary with spatial scale, because coexistence mechanisms are scale dependent. This scale dependence may alter the shape of the species-area relationship (SAR), which can be described by changes in species richness (S) as a power function of the sample area (A): S = cAz , where c and z are constants. We analysed the effects of experimental manipulations of nutrient supply and herbivore density on species richness across a range of scales (0.01-75 m2 ) at 30 grasslands in 10 countries. We found that nutrient addition reduced the number of species that could co-occur locally, indicated by the SAR intercepts (log c), but did not affect the SAR slopes (z). As a result, proportional species loss due to nutrient enrichment was largely unchanged across sampling scales, whereas total species loss increased over threefold across our range of sampling scales.

Negative effects of nitrogen override positive effects of phosphorus on grassland legumes worldwide.

Anthropogenic nutrient enrichment is driving global biodiversity decline and modifying ecosystem functions. Theory suggests that plant functional types that fix atmospheric nitrogen have a competitive advantage in nitrogen-poor soils, but lose this advantage with increasing nitrogen supply. By contrast, the addition of phosphorus, potassium, and other nutrients may benefit such species in low-nutrient environments by enhancing their nitrogen-fixing capacity. We present a global-scale experiment confirming these predictions for nitrogen-fixing legumes (Fabaceae) across 45 grasslands on six continents. Nitrogen addition reduced legume cover, richness, and biomass, particularly in nitrogen-poor soils, while cover of non-nitrogen-fixing plants increased. The addition of phosphorous, potassium, and other nutrients enhanced legume abundance, but did not mitigate the negative effects of nitrogen addition. Increasing nitrogen supply thus has the potential to decrease the diversity and abundance of grassland legumes worldwide regardless of the availability of other nutrients, with consequences for biodiversity, food webs, ecosystem resilience, and genetic improvement of protein-rich agricultural plant species.

Preference, performance, and chemical defense in an endangered butterfly using novel and ancestral host plants.

Adoption of novel host plants by herbivorous insects can require new adaptations and may entail loss of adaptation to ancestral hosts. We examined relationships between an endangered subspecies of the butterfly Euphydryas editha (Taylor's checkerspot) and three host plant species. Two of the hosts (Castilleja hispida, Castilleja levisecta) were used ancestrally while the other, Plantago lanceolata, is exotic and was adopted more recently. We measured oviposition preference, neonate preference, larval growth, and secondary chemical uptake on all three hosts. Adult females readily laid eggs on all hosts but favored Plantago and tended to avoid C. levisecta. Oviposition preference changed over time. Neonates had no preference among host species, but consistently chose bracts over leaves within both Castilleja species. Larvae developed successfully on all species and grew to similar size on all of them unless they ate only Castilleja leaves (rather than bracts) which limited their growth. Diet strongly influenced secondary chemical uptake by larvae. Larvae that ate Plantago or C. hispida leaves contained the highest concentrations of iridoid glycosides, and iridoid glycoside composition varied with host species and tissue type. Despite having largely switched to a novel exotic host and generally performing better on it, this population has retained breadth in preference and ability to use other hosts.

Opposing community assembly patterns for dominant and nondominant plant species in herbaceous ecosystems globally.

Biotic and abiotic factors interact with dominant plants-the locally most frequent or with the largest coverage-and nondominant plants differently, partially because dominant plants modify the environment where nondominant plants grow. For instance, if dominant plants compete strongly, they will deplete most resources, forcing nondominant plants into a narrower niche space. Conversely, if dominant plants are constrained by the environment, they might not exhaust available resources but instead may ameliorate environmental stressors that usually limit nondominants. Hence, the nature of interactions among nondominant species could be modified by dominant species. Furthermore, these differences could translate into a disparity in the phylogenetic relatedness among dominants compared to the relatedness among nondominants. By estimating phylogenetic dispersion in 78 grasslands across five continents, we found that dominant species were clustered (e.g., co-dominant grasses), suggesting dominant species are likely organized by environmental filtering, and that nondominant species were either randomly assembled or overdispersed. Traits showed similar trends for those sites (<50%) with sufficient trait data. Furthermore, several lineages scattered in the phylogeny had more nondominant species than expected at random, suggesting that traits common in nondominants are phylogenetically conserved and have evolved multiple times. We also explored environmental drivers of the dominant/nondominant disparity. We found different assembly patterns for dominants and nondominants, consistent with asymmetries in assembly mechanisms. Among the different postulated mechanisms, our results suggest two complementary hypotheses seldom explored: (1) Nondominant species include lineages adapted to thrive in the environment generated by dominant species. (2) Even when dominant species reduce resources to nondominant ones, dominant species could have a stronger positive effect on some nondominants by ameliorating environmental stressors affecting them, than by depleting resources and increasing the environmental stress to those nondominants. These results show that the dominant/nondominant asymmetry has ecological and evolutionary consequences fundamental to understand plant communities.

General destabilizing effects of eutrophication on grassland productivity at multiple spatial scales.

Eutrophication is a widespread environmental change that usually reduces the stabilizing effect of plant diversity on productivity in local communities. Whether this effect is scale dependent remains to be elucidated. Here, we determine the relationship between plant diversity and temporal stability of productivity for 243 plant communities from 42 grasslands across the globe and quantify the effect of chronic fertilization on these relationships. Unfertilized local communities with more plant species exhibit greater asynchronous dynamics among species in response to natural environmental fluctuations, resulting in greater local stability (alpha stability). Moreover, neighborhood communities that have greater spatial variation in plant species composition within sites (higher beta diversity) have greater spatial asynchrony of productivity among communities, resulting in greater stability at the larger scale (gamma stability). Importantly, fertilization consistently weakens the contribution of plant diversity to both of these stabilizing mechanisms, thus diminishing the positive effect of biodiversity on stability at differing spatial scales. Our findings suggest that preserving grassland functional stability requires conservation of plant diversity within and among ecological communities.